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ABSTRACT 

We perform a detailed study of the Globular Cluster (GC) system in the 
galaxy NGC 5866 based on F435W, F555W, and F625W (~ B, V, and R) HST 
Advanced Camera for Surveys images. 

Adopting color, size and shape selection criteria, the final list of GC can- 
didates comprises 109 objects, with small estimated contamination from back- 
ground galaxies, and foreground stars. 

The color distribution of the final GC sample has a bimodal form. Adopting 
color to metallicity transformations derived from the Teramo-SPoT simple stellar 
population model, we estimate a metallicity [Fe/H]~ —1.5, and —0.6 dex for the 
blue and red peaks, respectively. A similar result is found if the empirical color- 
metallicity relations derived from Galactic GCs data are used. 

The two subpopulations show some of the features commonly observed in the 
GC system of other galaxies, like a "blue tilt", higher central concentrations of 
the red subsystem, and larger half-light radii at larger galactocentric distances. 
However, we do not find evidence of a substantial difference between the average 
sizes of red and blue clusters. 

Our analysis of the GC Luminosity Function indicates a V-band Turn-Over 
Magnitude V^ OA/ =23.46±0.06, or M™ M ~ -7.29 ±0.10 mag, using the distance 
modulus derived from the average of SBF and the PNLF distances. The absolute 
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Turn-Over Magnitude obtained agrees well with calibrations from literature. The 
specific frequency is measured to be Sn = 1.4 ± 0.3, typical for galaxies of this 
type. 

Subject headings: globular clusters: general — galaxies: individual, NGC5866 
— galaxies: star clusters — galaxies: distances 



Introduction 



The study of the Globular Cluster (GC) system in galaxies is one of the fundamental 
keystones to understand the physical processes at the base of the formation and evolution of 
galaxies. The investigation of GC system properties involves aspects of several different fields 
in astronomy: stellar evolution, hierarchical assembly, the determination of extra-galactic 
distances, etc. 

The first applications of the GCs as an investigation tool came in the earliest years of 
extra-galacti c astronomy . Single bright GC were recognized as possible distance indicator 
in the 1950s (jBaumlll955l ). but it took few decades to develop the technique as it is applied 
now. The systematic study of the GC Luminosity Function (GCLF), in fact, has shown 
that it has a universal nearly gaussian shape. It is now recognized that the Turn Over 
Magnitude (TOM) of the GCLF is a reliable distance indicator for most galax ies, although 
some exceptions exist (IFerrarese et al.ll2000l ; lRichtlerll2003l ; I Jordan et al.l 120071 ). 



Another, more recent discovery in the field of GC astronomy, is that the presence of a 
color bimodality is quite common in galaxies. The color bimodality has been mos tly inter- 
preted as a bimodal metallicity distribution (IHarrid l200ll : iBrodie fc Straderll2006l ). leading 
to develop several possibl e GC formation scenar ios, useful to constrain the history of galaxy 
formation. More recently, iKundu &: Zepj (120071 ) using near-IR photometry, have shown that 
the GC bimodality in the giant ellipticals M87 and NGC4472 might represent true metallic- 
ity subpopulations. However, whether or not the bimodal color distribution is linked to the 
presence of two metallicity subpopulations is still a matter of debate, and there are claims 
that the optical properties of a bimodal GC system are not necessarily due to a bimodal 
metallicity distribution. For example, it has been shown that a bimodal color distributions 
might be generated by the combination of a unimodal metallicity distribution and a no n- 
linear color- metallicity relation of the appropriate shape ( lRichtlerll2006l ; lYoon et al.ll2006l ). 



Furthermore, the advent of space-based astronomy has provided new information in this 
field. Most of all, the high resolution optical imaging of HST makes it possible to examine the 
shape properties, i.e. size and ellipticity, for an increasing number of extra-galactic GCs, up 



- 3- 



to distances as far as the Virgo Cluster (1 Jordan et al.l 120071 ) and beyond. The examination 
of the characteristic sizes, and radial distribution of the blue and red GC subpopulations, 
enables the study of new phenomen a related to the formation history of the GC system 
f lLarsen fc Brodie!l2003l : IjordaiJliooi l. 



In this paper we take advantage of archival HST images of the galaxy NGC5866, taken 
with the Advanced Camera for Surveys (ACS) in the Wide Filed Camera mode, to carry out 
the first detailed multi-wavelength study of the GC system in this galaxy. 

NGC5866 is a lenticular galaxy, with a sharp dust lane nearly edge on. Being a disk 
galaxy, the GC system is expected to be less popul ated with respect to ellipticals (e.g. 



Harris &: van den Bergh] Il98ll ; iBrodie fc Strader 



20061 ). However, the quality of these ACS 



data allows us to analyze in some detail the clusters properties. 



The weighted average of distances moduli obtained usi ng the SBF (|Tonrv et al.l 12001 
using the new zeropoint from Jensen et al, 2003), and PNLF (jCiardullo et al.ll2002l ) methods, 
provide a distance modulus /i = 30.75 ± 0.08, or 14.1±0.5 Mpc. At this distance, the GCs 
of the galaxy are near the ACS resolution limit, allowing us to carry out the size and shape 
analysis of the GC system. 

The analysis of GC properties will be presented in the next sections. The observations 
and data reduction procedure are introduced in Section 2. Section 3 deals with the GC 
photometry and selection criteria. The color bimodality, color-magnitude relations, and 
the spatial distribution of the total GC system, as well as the two GC subpopulations, are 
discussed in Section 4. In the same Section we also make a comparison with stellar synthesis 
models. Finally, a study of GCLF and specific frequency is also presented in Section 4. We 
summarize and conclude the paper in Section 5. 



2. Data reduction 



The ACS/WFC camera observations of NGC5866 were taken in February 2006 with 
the F435W, F555W and F625W filters (~ B, V, and R). The total exposure times for the 
three different passbands were: 3900s, 2800s and 2200s, respectively. 

Stan dard reduction was c arried out on raw images using the APSIS image processing 
software (IBlakeslee et al.ll2003l ). After cosmic-ray rejection, alignment, "drizzling", and final 
image combination, we adopted an iterative procedure to determine the sky background and 
the galaxy best fit models, independently for the three passbands. 



Since the galaxy fills the whole field of view of the ACS, the sky background determined 
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in the outer corners of the CCD image is not accurate. We obtain the sky values extrapolating 
from the fit of a de Vaucouleurs r 1 / 4 profile. 

Briefly, after masking out a few bright/saturated foreground stars, we adopt a pro- 
visional sky value from the corner with the lowest median number of counts. This value 
is subtracted to the image before starting the galaxy fitting procedure. Then, we fit the 
galaxy isoph qtes using the IRAF /STSDAS task ELLIPSES, which is based on the method 
described by iJedrzejewskil (119871 ). Once the preliminary galaxy model is subtracted from 
the sky-subtracted frame, a wealth of faint sources appears together with the dust disk. A 
mask of these sources is obtained using a sigma detection method, and combined with the 
previous mask. The new mask is then fed to ELLIPSE to refit the galaxy's isophotes. 

After the geometric profile of the isophotes has been determined, we measure the galaxy 
flux in elliptical isophotes. At last, to improve the estimation of the sky, we fit the surface 
brightness profiles with a de Vaucouleurs r 1//4 profile plus a constant sky offset. The sky offset 
fitted is then adopted as the new sky value, and the whole procedure of galaxy fitting, source 
masking, surface brightness profile fitting and sky estimation is repeated, until convergence. 
This procedure is adopted for all three filters. 

Once the sky and galaxy model have been subtracted from the frame, some large-scale 
deviations are still present in the frame. We remove these deviations using the background 
map obtained running SExtractor (IBertin fc Arnoutslll996l ) on the sky and galaxy subtracted 
frame. Part of the procedure adopted here is the same developed by ICantiello et al.l (120051 ) 
to process ACS images for SBF analysis. We refer the reader to that paper for more details 
on the whole procedure of sky+galaxy+large scale residuals subtraction. 



It is worth noting that in ICantiello et al.l (120051 ) we have made some checks on the 
stability of the sky estimation against the fitting profile chosen. There, it was shown that, 
adopting a more general Sersic r 1//n profile, the sky values obtained agree wi thin uncertain- 



ties w ith those obtained adopting a de Vaucouleurs profile. In addition, in ICantiello et al. 



(|20070. we compared the sky values obtained applying our method with the sky obtained by 
Sikkema et al.l (120061 ) using a different approach, over a sample of 6 peculiar galaxies. As a 
result we found that the differences between the two sky background determinations is on 
average ^3%. 



1 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the As- 
sociation of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 
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Figure [T] shows the I-band ACS image of NGC5866_| (left panel), and a zoom of the 
sky+galaxy+large scale residuals subtracted frame (right panel, "residual" frame hereafter), 
the nearly edge-on dust disk in the galaxy is easily recognizable in this second panel. 



3. Globular Cluster selection and photometry 



The residual images were used for the detection and photometry of point-like and ex- 
tended sources present in the frames. The source photometry was obtained independently on 
all three frames using SExtractor. After source detection, the output catalogs were matched 
using a 0.1" radius. 



To derive the aperture correction, standard growth curve analysis (jStetsonlll990l) was 
perfo rmed on selected, isolated compact sources. Extinction corrections (values from lSchlegel et al 
19981 ) and aperture corrections were applied before transformin g the ACS magnitude s into 
the standard B, V and R magnitudes, using the prescriptions by lSirianni et al.l (120051 ). The 
transformations from ACS to sta ndard BVRI sy stem magnitudes were obtained adopting 
the synthetic transformations from ISirianni et al.l We have checked the compat ibility of these 



magnitudes with the ones derived adopting the empirical transformations from lSirianni et al. 



As a result, we have found that the differences between the magnitudes obtained with the 
two set of transformations agree on average within 0.02 mag, at least in the range of colors 
of GC candidates which will be adopted later on. 



We have chosen to select GC candidates using an approach similar to ISpitler et al. 



( 120061 . S06 hereafter), based on the color, shape and size of the objects. At the distance of 
NGC 5866, the typical GC half-light radius, ~3 pc, is at the limit of the angular res olution 
of th e ACS camera (0.05 arcsec/pixel). Under such condition, the software ishape ( jLarsen 
19991 ) can be used to estimate the intrinsic shape parameters of slightly resolved objects. To 
run the ishape software an oversampled PSF of the original image is needed. Based on their 
colors, we have selected a sample of star candidates from each frame and derived the PSF 
using the photometry package DAOPHOT, in the IRAF environment. The objects selected 
for PSF analysis were a posteriori confirmed to be stars. 

In order to save computational time, we have removed from the initial catalog of matched 
sources all the sources easily recognizable as background galaxies. The new catalog was then 
fed to ishape, with the PSFs, to derive the shape parameters of the listed objects in the three 
different passbands. One of the main parameters to run ishape is the convolution kernel, 



2 Also visit the web site 



http://heritage.stsci.edu/2006/24/index.htmlfor a nice color picture of the galaxy. 
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which, once convolved with the PSF, is used to fit the ellipticity, orientation, FWHM, etc. 
of the object. After various tests, we found that the best results were obtained adopting the 
analytic King profile, with concentration parameter c=30. 

The final shape parameters for each object were then adopted as the average between 
the three estimations from the B, V and R frames. Those objects with half-light radius in 
one band more than 3-cr different from the V-band value, were rejected from the catalog. 
Similarly to S06, we have repeated such rejection procedure twice on the whole catalog of 
objects. 

Next, we have rejected those object with ishape ellipticities > 0.5 (< 10 objects at 
all). As a final selection criteria we have used the (B-R) and (B-V)o color of the remaining 
sources. Only those objects whose color fell within the intervals 0.9< (B — R) < 1.7, and 
0.5 < (B — V) < 1.1 were taken as globular cluster candidates. After these selections, the 
final catalog of GC candidates consisted of 109 objects. The list of GC candidates, and their 
photometric, spatial and shape properties are reported in Table [H 

Before using these data, it is useful to estimate the amount of contamination from 
background galaxies and foreground stars. 

Concerning the contamination due to background galaxies, given the similarity of NGC 5866 
data with those of the Sombrero galaxy by S06, we adopt the conclusions from these au- 
thors. S06 have reduced and analyzed B and V images from Hubble Deep Field-North of 
the GOODS survey, and found that the number of background objects matching the GC 
selection criteria is ~ 4% the final Sombrero GC sample. When this number is scaled to the 
imaging area of NGC 5866 (~l/6 of the Sombrero), and to the number of objects classified 
as GC, we find that in our case the contamination due to background galaxies is ~ 6% the 
final number of GCs. 

The number of galactic foreground st ars expected in the color ranges adopted for GC 



selection have been calculated using the iBahcall &: Soneiral (Il980l ) model. A number of 



approximately ~ 8 stars is predicted to fall in the range of colors adopted for the GC 
selec tion. We fin d a sim ilar, but somewhat lower number (~ 5-6 stars) using the models 



from iRobin et al.l (120031 ). These numbers have to be compared with the number of objects 
actually matching with the adopted color ranges, and rejected because of the small half-light 
radius, that is 4 objects at all. Given the small number statistics, we consider the number of 
contaminating stars from models consistent with the observed one, thus we can safely assume 
that the final sample of GC in Table [T] is likely free of contamination from foreground starg^l 



3 If the expected and the observed number of stars are taken as they are, i.e. without taking into account 
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Discussion 



4.1. Color histograms 

Figure [2] shows the (B-R) , (B-V) and (V-R)o color histograms for the 109 GC candi- 
dates. To_Je^t_the_j2imodah^y of the GC sample we use the KMM 
code (jMcLachlan fc Basfordlll988l ; lAshman et al.lll994l ). This code evaluates the likelihood 
whether or not the color distribution of sample is better represented by two gaussians or one 
single gaussian, the positions of the peaks are evaluated as well. 

Running the KMM code on the (B-R) color, returns a likelihood of ~96% that the GC 
color distribution is bimodal, with peaks at (B-R) =1.12 (blue peak), and (B-R) =1.35 (red 
peak). The GCs are evenly assigned to each group. 

An independent KMM run on the (B-V)o color, returned a slightly lower but significant 
bimodality probability (93 %), and peaks at (B-V)o=0.67, (B-V)o=0.84 mag. 

The narrow range of GCs (V-R)o color, together with the small separation between the 
two possible peaks (< 0.1 mag for the ~ 600 GCs of the Sombrero galaxy), and the relatively 
small number of GC candidates, hampers any meaningful solution by KMM for this color. 

At pres ent, it is still a matter of debate what is the origin of the bimodality. The recent 



review from iBrodie fc Straderl (120061 ) listed some of the different scenario proposed (major 
mergers, dissipationless accretion, hierarchical merging, etc.). In general the photometric 
data available, and the few spectroscopic measurements seem to point out that the primary 
difference between the two GC subpopulations is due to metallicity, without the need of 
invoking age variations. On the other hand, recently, it has been shown that a nonlinear color- 



Richtler 


2006; 


Yoon et al. 


2006) 



If we assume the picture where the two color peaks are generated by a bimodal metal- 
licity distribution, the metallicities corresponding to the peaks can be derived by assuming 
so me [Fe/H] versus color r elations. Adopting the view that GCs are mostly old, and using 
the Raimondo et al. ( 2005 . R05 hereafteiB) models in the metallicity range —1.8 < [Fe/H] 
< 0.3, for ages in the interval 9 < t (Gyr) < 14, we derive [Fe/H] = (3.67 ± 0.13) (B - 
R)o — (5.60 ± 0.18). We have verified the reliability of this equation by comparing the the- 



the uncertainties, we are lead to conclude that ~ 3% of the final GC sample is composed by contaminating 
stars. 



lr The R05 models are available at the Teramo-SPoT group website: www.oa-teramo.inaf.it/SPoT 
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oretical [Fe/H]^. values predicted for the Galactic GCs, with the observational [Fe/H] {, s . 
estimations. For t his com p arison we have adopted the (B-R) and [Fe/H] obs values from 
the updated online Harris! ( 19961 ) catalog of Galactic GCg§. As a result we find that the 



median difference is ([Fe/H] t h. — [Fe/H] b s .) = —0.02 ± 0.16 dex, confirming the reliability 
of the adopted color-metallicity relation. The good agreement of the models with Galactic 
GC data can also be seen in Figure [3] (left panel) where we show the fit of the [Fe/H] versus 
(B-R)q from the R05 models, and the Galactic GCs data. 

Using the above [Fe/H] versus (B-R) equation, the NGC5866 peak metallicities are 
[Fe/H]=— 1.48 ± 0.23, and — 0.64±0.25 dex for the blue and red GC peaks, respectively. 

We have repeated the above steps using (B-V)o colors instead of the (B-R) . The color- 
metallicity relation in this case is [Fe/H] = (5.48±0.19)(B-V) - (5.13±0.16) (Fig. El right 
panel), which implies [Fe/H]=— 1.49±0.20, and — 0.54±0.23 dex for the two subpopulationg^l 

These metallicity values can be considered quite average in bimodal GC s y stems . In 



fact, adopting the estimations for other galaxies reported by iBrodie fc Straderl (120061 ). we 
find on average [Fe/H]=— 1.47±0.12 and — 0.48±0.16 dex for the sample of normal elliptical 
and SO galaxies. 



4.2. Blue tilt 



Inspecting the color magnitude diagram of GCs, it has been shown that in som e giant 
galaxies a color- luminosity correlation exists for the blue GCs ([Harris et al.l 120061 ) . The 
interpretation of this correlation must be taken carefully, due to the limits of the GC selection 



criteria. Moreover, this phenomenon has not been detected in all galaxies (IMieske et al.ll2006 
Brodie fc Straderl l2006h 



To explain the existence of the color-luminosity correlation two different GC formation 
scenarios have been proposed: the self-enrichment, and the pre-enrichment one. The dif- 
ferences between these two scenarios is that in the first one the proto-GCs are able to trap 
and absorb the metals formed by stars within the GC; in the pre-enrichment scenario it is 
assumed that higher mass proto-GCs come from intrinsically higher metallicity clouds - see 
Strader et al.l (120061 ) and references therein for more details. 



'The catalog can be found at the web address: http://www. physics. mcmaster.ca/^harris/mwgc.dat 



6 Using the empirical color-metallicity relations derived by S06 from Galactic GC data, the GC color peaks 
correspond to [Fe/H]~ —1.5, and [Fe/H]~ —0.75 dex for the blue and red subpopulations respectively. 
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Figure H] shows the color magnitude diagrams for the NGC5866 GCs. For the (B-R)o 
and (B-V)o we independently used the KMM output classification to separate the GC system 
into subpopulations. For the (V-R) we used the combined results from the above colors to 
separate between the two subpopulations. It is worth noting that the sharp separation limits 
adopted to discriminate between the two sub-populations introduce a bias in the results, since 
real GC subsystems must overlap in some color range. Keeping in mind this warning, in 
Figure H] we show the linear least-squares fits to the single subsystems. 

We find a slight color-luminosity correlation for the blue GCs: the slope of the (B-R)o 
color versus R magnitude relation is —0.012 ± 0.009. The blue tilt is more evident if (B-V)o 
and (V-R)o data are taken into account^, since with slopes — 0.017±0.009, and — 0.013±0.003, 
respectively. The slopes of the color-luminosity relations of the red sub-population are in all 
cases consistent with zero. 

Using the color metallicity relations reported in the previous section to convert the tilt 
in terms of metallicity differences, we find that the brightest blue GC have on average ~ 0.4 
dex higher [Fe/H] respect to the faintest ones. 



It is interesting to note that the study of a sample of bright ellipticals lead lHarris et al. 



( 120061 ) to conclude that color-luminosity relation appears more clearly in galax ies with larger 



number of clusters. An analogous result has been found by iMieske et al.l (120061 ) by analyzing 
the GC properties for a sample of 79 Virgo Cluster galaxies. A similar conclusion can be 
drawn by comparing our data with the results of the Sombrero galaxy by S06. In fact, 
NGC 5866 is morphologically similar to the brighter Sombrero galaxy, but its GC population 
within ~1.5 galactic effective radii is ~ 1/5 the GC population in ~1.5 effective radii of the 
Sombrero^, and the (B-R) versus R slope for blue GCs is almost one third the one derived 
from the Sombrero galaxy. 



4.3. Spatial properties 

One of the common results concerning the projected radial density of the GC subpopula- 
tions in galaxies, is that they can be fitted by a power law, with the red subpopulation being 
more centrally concentrated than the blue one. More specifically, it has been shown that the 



7 For this Figure we have determined the blue and red peaks of the (V-R)o color using: (V-R)o,biue = (B- 
R.)o,fei«e-(B-V)o,6i« e =0.45 mag, and (V-R)o,red=(B-R)o,red-(B-V)o,red=0.51. These numbers agree with the 
observed (V-R)o peaks in the Sombrero (S06). 

8 NGC 5866 images have nearly the same absolute magnitude completeness limit of the Sombrero images. 
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outer parts of GC systems are typically well- fit by power laws, while the central few kpc. 



hidde n in these NGC5866 data, often show the presence of cores (see, e.g., lAshman &: Zepi 



19981 ). 



Figure [5] shows the GC surface density in NGC 5866 versus the logarithm of galacto- 
centric distance, \ogR gc . The upper panel shows the total density profile, while in the lower 
panel the profiles of the two GC subpopulations are shown separately. From this figure it 
is readily recognizable that i) the red subpopulation appears slightly more centrally concen- 
trated respect to the blue one, ii) the density profiles are well fitted by power laws, with 
slopes: —1.7 ± 0.3, —0.9 ± 0.3, and —1.3 ± 0.2 for the red, blue, and total populations re- 
spectively. By comparing the slope of the total density profile with data from literature we 



find th at NGC 5866 has an average behavior. In fact, using the compilation by iKissler-Patig 



( 119971 ). we find that slopes in the range [—2.5,-1] are quite common in galaxies of total 



magnitude similar to NGC 5866, i.e., Mst ~ —20 mag. 

The radial profiles of the (B-R)q color for the two subpopulations do not show significant 
radial gradients (upper panels in Figure E]). However, due to the different concentrations of 
the red and blue GCs shown in Figure El the total cluster population shows a clear color 
gradient (Figure [6l lower left panel). In the lower right panel of Figure [6] we show the (B-R)o 
color profile of the galaxy, the best fit line to the total GC color profile is shown, too. As 
can be recognized from this figure, the color profile of field stars closely resembles the color 
profile of the total GC system, with a ~ +0.07 mag shift, i.e. ~ 0.3 dex higher metallicity. 

These observational properties can be interpreted with a scenario where the metal poor 
blue GCs form in the early stages of the galaxy life, while the red GC component and the 
field stars are formed later on in a metal enriched potential well. Or, it could just be a 
metallicity gradient for both the GCs and field stars, combined with the nonl inear color- 



met a l ii city relation that originates the observed GC bimodal color distributions (jYoon et al. 



20061 : lRichtledl2006h 



One other issue in understanding GC system formation and evolution, is the relation 
existing between the average sizes of GC subpopulations, and the variation of half-light radii 
with galactocentric distances. 

Some observational data have shown that the half-light radii of the red subpopulation are 
on average 10% to 20% smaller than blue G C. Two different explan ations have been given 



to this phenomenon. The one proposed by lLarsen fc Brodiei (120031 ) involves a projection 



effect, due to the different spatial distribution of red and blue GCs, together with a strong 
correlation of the GC size with the galaxy radius. This model predicts that the size differences 
should be larger in the inner galactic regions, and disappear in the outskirts of the galaxy. In 



the other scenario, proposed by I Jordan! (120041 ). mass segregation and stellar evolution effects 
are the leading causes of the blue/red GC size differences. In this scenario little change of 
relative GC sizes with galactocentric distances is expected. 



By comparing the average half-light radii of the red and blue GC in NGC 5866 we do not 
find any significant difference between the two subpopulations. Furthermore, the size versus 
galactocentric radius comparison, shown in Figure [7J reveals a slight tendency of outer GC 
to have larger half-light radii. 

Due to the small number of GCs available, the absence of a clear size difference between 
the red and blue GCs, and the uncertainties in the half-light radius versus R gc correlation, 
these observations cannot be used to clearly support either the mass segregation or the 
projection scenarios. 



4.4. Comparison with models 

Figure [8] shows a comparison between the observed colors of NGC 5866 GCs and the 
predictions from R05. Each panel in the figure shows a Simple Stellar Population (SSP) 
model with a defined [Fe/H], and for ages 1 < t (Gyr) < 14. 

At first glance, these panels will give little information about the age and metallicity 
properties of the observed GCs, due to the strong overlap of models with different phys- 
ical properties. This is not related to the particular choice of SSP models. In order to 
reduce model systematics , we have also considered a data to models comparison using the 



Bruzual fc Charlotl (120031 ) SSP models. The new set of models essentially agrees with the 
R05 one, and the overlap between models at different [Fe/H] is still present. Moreover, since 
the R05 models are optimized to match in detail the observational features of Color Mag- 
nitude Diagrams of Galactic and Magellanic Clouds stellar clusters, in this section we will 
take the set of stellar synthesis models from R05 as reference ones. 

If no constraint is put on models, the conclusions that can be drawn from the comparison 
shown in Figure [8] are mainly: i) SSP models, i.e. single age and single metallicity stellar 
systems, reproduce nicely the observational color-color properties of GCs; ii) except one 
single object, no old super metal-rich (t > 7 Gyr, [Fe/H]>[Fe/H] ) GC seems to be present; 
Hi) the reddest GCs in the sample [(B-R) ~1.6, (B-V) ~1] are either old, t~ 14 Gyr, solar 
metallicity objects, or younger, t~ 5 Gyr, super-solar metallicity stellar systems. 

However, if we limit the models to the range of old age models (9 < t (Gyr) < 14) some 
stronger constraints on the age and metallicity of the GC can be obtained. To support this 
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choice, as mentioned before, we recall that old GC ages have been found in most of the GC 



systems in galaxies both with photometric and spectroscopic data (IBrodie & Straderll2006l ). 
In addition to this, as shown in Figure M (upper panels), the range of (B-R) and (B-V) 
colors for NGC 5866 and for Galactic GCs appears to be very similar. If we assume that such 
color similarity is also associated with common physical characteristics of the GC systems, 
than we can adopt for NGC 5866 an age rang e similar to the one of Galactic GC, that is 9 



< t (Gyr) < 14 (e.g. iDe Angeli et al.ll2005h . 



Based on these considerations, we can go back to the comparisons shown in Figure El 
taking into account only the models in the selected age interval. With this constraint we 
find that the blue tail of GCs has [Fe/H]< —1.8 dex. If we extrapolate the [Fe/H] versus 
(B-R)q equation derived from R05 models using the same age limitations adopted here, then 
the most metal poor GC in NGC 5866 has [Fe/H]~ —2.3 dex. On the other hand, inspecting 
the color properties of red clusters, the conclusion is that the reddest GCs in NGC 5866 are 
old objects, t~14 Gyr, with nearly solar metallicity. 

Before concluding this section, we find interesting to note that if the same arguments 
are applied to the GC system in the Sombrero galaxy (data from S06, shown in the right 
panels of Figure Mj, the metallicity range spanned by the GC system is —2.2 ^ [Fe/H] ^ 
0.6. Thus, while the blue tails of the GC system appear to be quite similar in the three 
galaxies, the red tail of the GCs in the Sombrero includes also a population of very metal 
rich clusters. 



4.5. Globular Cluster Luminosity Function 



The review by iRichtlerl (120031 ) gives a fairly detailed study on the use of GCs as dis- 
tance in dicators. By comparing GCLF distances with those derived using the SBF method, 
Richtlerl proves that the TOM is a good distance indicator, although some exceptions exists, 
related to a possible contamination from intermediate age GCs, or young massive clusters. 

Figure [10] shows the V-band completeness functions for NGC 5866 derived within three 



annular regions with radii R, 



■gc 



arcsec 



< 35, 35 < Rg C (arcsec) 



< 70, and R. 



',)<- 



arcsec] 



> 70. 



A 95% completeness level is reached at Vo = 23.9, 24.8, 25.1 mag for the three different 
annuli, respectively. The completeness functions have been derived via standard artificial 
star tests, using the DAOPHOT task addstar to generate the artificial stars, then we ran 
SExtractor for the detection of sources, as described in Section 2. 



Adopting SBF distances. IRichtlerl (120031 ) obtained an absolute TOM M Vfi = -7.35±0.25 
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ma^|, which at the distance of NGC5866 translates into an apparent TOM V™ A/ ~ 23.46 
mag. This means that the GCLF is fairly complete up to ~0.5, 1.4, and 1.7 magnitudes 
fainter than the TOM in the three different annular regions taken into account. 

Fitting a gaussian distribution to the GCLF corrected for incompleteness, gives a TOM 
V™ A/ = 23.46 ± 0.06 mag, with a dispersion a = 1.13 ± 0.05. Adopting the SBF+PNLF 
average distance /i = 30.75 ±0.08, th e absolut e TOM is M^ M = -7.29 ±0.10 mag, in nice 



agreement with the value provided by iRichtler . 



4.6. Specific Frequency 

Although the field of view of the ACS camera does not sample the whole NGC5866 
area, we have obtained an estimation for the Specific frequency, Sn, of GCs over the ACS 
sampling area. 



The concept of specific frequency, Sn = A' r Gc , xl0 ' 4< - Ml/+15 ' ) , was introduced by lHarris fc van den Bergh 



( I198ll ) as a tool to quantify the amount of stars in GCs with respect to the amount of field 
stars in a galaxy. To derive Sn we choose to not do any extrapolation to derive the total Nqc 
and the total magnitude, m v , to avoid the un certainty introduced b y the profile extrapola- 
tion. Instead, adopting an approach similar to lBlakeslee et al.l (119971 ). we adopt the number 
of GCs found in the ACS sampling area, after completeness correction, and the total mag- 
nitude obtained from the same masked frame that we used to detect GCs [my = 10.9 ±0.1 
mag). 



As noted by lBlakeslee et al.l (119971 ). if the GC system follows the same radial profile of 
the galaxy light in the outer regions, there should be little diff erence between t he global Sn 



and this Sn over a limited area, referred to as "metric Sn" by lBlakeslee et al 



With these assumptions, we derive Sn = 1.4 ±0.3, which is wit hin the normal Sn- range 
for galaxies in the same luminosity class of NGC5866 ( S m =0.7-7.5. lKissler-Patiglll997l ). We 
have also estimated the total Sn extrapolating the total number of GCs. By integrating 
over the surface density, corrected for magnitude completeness, we obtained Nq^"- 1 ~ 400. 
Adopting a total magnitude m v ~ 9.8 ± 0.1 mag from the RC3 catalog, we find Sj ot ' n - 
1.6 ± 0.4, in nice agreement with the value found over the limited ACS area. 



'N 



Finally, since the GCs are evenly distributed between the two subpopulations, we find 



9 We have corrected th e TOM value reported by Richtler ( 20031 ) using the SBF zeropoint correction as 
discussed bv I Jensen et al.l f|2003h . 
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that the specific frequency for the blue and red subsystems is S^ ue ~ S r ^ d ~ 0.7. 

5. Conclusions 

We have used archival HST data to study the GC system of NGC 5866 galaxy in the B, V 
and R passbands. Using standard color, size and shape selection criteria, the final list of GC 
candidates consists of 109 objects. The estimated contamination from background galaxies 
is ~ 6 % of the sample, while essentially no contamination is expected from foreground stars. 

The color distribution of the GC system has a bimodal form. Adopting the color- 
metallicity relations from R05 models, we find that the blue and red peaks correspond to 
a metallicity [Fe/H]~ —1.5, and —0.6 dex, respectively. These metallicity peaks are quite 
average among the galaxies with known color bimodality. Similarly, we have found that the 
specific frequency of GCs is Sn = 1-4 ± 0.3, which is a normal value for galaxies of the same 
luminosity class of NGC 5866. 

By inspecting the color magnitude diagram of the two GC subsystems, a color-luminosity 
correlation appears in the blue GCs, with brighter objects being slightly redder than the 
fainter ones. Interpreting this feature as a luminosity-metallicity correlation, the bright blue 
GCs are found to be 0.4 dex more metal rich than the fainter ones. 

Inspecting the spatial properties of the GCs, we find that i) the red subpopulation 
appears to be more centrally concentrated than the blue one; ii) the radial density of the 
total GC population, and of the two subpopulations taken separately, follows a power law 
radial profile; and that Hi) the GCs have larger half-light radii at larger galactocentric 
distances. We do not find any sign of a significant difference between the average sizes of 
the subpopulations. 

Taken separately, the two GC subsystems do not show significant color gradients. How- 
ever, due to the different radial concentrations, the total GC system shows a color gradient, 
very similar to the one of field stars, except that the GCs are systematically ~0.1 mag bluer 
in (B-R)o- If the GC to field stars color offset is simply associated to a metallicity difference, 
then the population of field stars at a given galactocentric distance is, on average, 0.3 dex 
more metal rich of the GC population at the same radius. 

A comparison of data with the R05 models has shown that no old, very metal rich (t>9 
Gyr, [Fe/H]>0) cluster is present. Moreover, if only models in the age range 9< t (Gyr) < 
14 are taken into account, then the metallicity range for the GC system in this galaxy is 
found to be —2.3 < [Fe/H] < 0.0. When compared to Galactic and Sombrero data, we find 
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that NGC 5866 and Milky Way GC share similar metallicity range, while the Sombrero has 
a tail of very metal rich GCs ([Fe/H] ~ 0.6). 

The study of the GC luminosity function has shown that the GC sample is fairly com- 
plete out to ~1.2 mag below the expected TOM. A gaussian fit to the GCLF provided us 
yTOM _ 23.46±0.06 mag. Adopting a distance modulus /zo = 30.75±0.08 from the weighted 
average of SBF and PNLF dis t ances , implies M^ M = -7.29 ± 0.10 mag, in agreement with 



the mean value from iRichtlerl (120031 ) . 



We thank the anonymous referee for helping us to improve this paper with useful sug- 
gestions. This work was supported by the NASA grant AR-10642, and by COFIN 2004 
under the scientific project "Stellar Evolution" (P.I.: Massimo Capaccioli). 
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Fig. 1. — Left panel: NGC5866 galaxy image from the HST ACS. The rectangle shows the 
location of area zoomed in the right panel. Right panel: A zoom of the disk of dust in 
NGC 5866 as seen in the residual frame. 
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Fig. 2. — Color histograms of the final sample of GC candidates. 
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Fig. 3. — Left panel: Galactic GC (B-R) color and [Fe/H] data (full dots) compared with 
the linear fit to the R05 models (dotted line, see text). The fit to the models nicely agrees 
with observational data. Right panel: as left panel, but for (B-V) color. 
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Fig. 4. — Color-Magnitude Diagrams of the final sample of GC. There is some evidence of the 
color-magnitude tilt for the blue GC subpopulation (empty squares). No significant tilt is 
found in the red peak GC candidates (full squares). The dashed lines show the least-squares 
fit to the data. 
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Fig. 5. — Upper panel: The surface density of the entire GC population. The least-squares 
profile fit (solid line) shows that the radial density profile is well reproduced by a power law, 
log ogc — (—1-3 ± 0.1) log R gc + (1.0 ± 0.1). Lower panel: Same as upper panel but the two 
subpopulations are shown separately. The dotted (dashed) line represents a linear fit to the 
data of red (blue) GCs. The fitted lines show that the red subpopulation is more centrally 
concentrated than the blue one. 
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-23- 




Fig. 9. — Color-color comparisons (upper panels), and (B-R) color histograms (lower panels) 
for Galactic, NGC 5866 and Sombrero clusters, as labeled. The least-squares fit lines are also 
shown. The Sombrero GC system includes a rich component of red, super metal rich globular 
clusters which is practically missing in the other two galaxies system. On the other hand, 
the blue tails of the GC distributions are quite similar for the three galaxies. 



-24- 



i i i i i i i i i i i i i i i i i i i r 



n A A 



A 

□ A 
□ A 

□ 



n □ 

K n A 



K □ 1 



a u a 

«a a 
□ □ 

□ 



H A 



R ^35" 



□ 35" < R gc 570" 
a R gc > 70" 



: H °n a 



° A^ 



K K St Xt K &4 

J I I I I I I I I I I I I I I I I I I I I I I l_ 



23 



24 



25 26 



27 



28 



Fig. 10. — V-band completeness functions within three different galactocentric annuli. 
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Fig. 11. — The observed GCLF histogram (corrected for incompleteness) and the best fit 
gaussian line. The TOM obtained is V™ M = 23.46 ± 0.06 mag (dashed line), dispersion 
a = 1.13 ± 0.05. Adopting such TOM magnitude we derive the absolute TOM M^° M = 
— 7.29 ± 0.10 mag, in agreement with the "universal" value from iRichtlerl ( 2003 ). The 95% 
completentess magnitude limits at different galactic annuli are also shown in the figure with 
vertical arrows. 
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Table 1. Globular Clusters in NGC5866 
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40 





84 





56 


2 


02 


-0 


49 


± 


0.25 


15 06 


26 


96 


55 


46 


16 


29 


37 


94 


24 


32 


± 





05 


23 


58 


± 





04 


23 


03 


± 





03 


1 


29 





74 





55 


5 


73 


-0 


95 


± 


0.25 


15 06 


26 


02 


55 


46 


18 


00 


38 


51 


22 


95 


± 





03 


22 


35 


± 





02 


21 


92 


± 





02 


1 


04 





60 





44 


2 


60 


-1 


81 


± 


0.19 


15 06 


26 


51 


55 


45 


52 


59 


38 


87 


22 


68 


± 





03 


21 


87 


± 





02 


21 


36 


± 





02 


1 


31 





81 





51 


2 


50 


-0 


74 


± 


0.20 


15 06 


25 


81 


55 


45 


17 


43 


40 


51 


24 


24 


± 





03 


23 


59 


± 





03 


23 


15 


± 





02 


1 


09 





65 





45 


2 


63 


-1 


59 


± 


0.20 


15 06 


27 


12 


55 


46 


15 


90 


41 


83 


24 


16 


± 





04 


23 


28 


± 





03 


22 


76 


± 





03 


1 


40 





87 





52 


2 


21 


-0 


41 


± 


0.24 


15 06 


27 


78 


55 


46 


26 


76 


41 


89 


23 


07 


± 





03 


22 


12 


± 





02 


21 


59 


± 





02 


1 


48 





95 





53 


2 


82 


-0 


03 


± 


0.21 


15 06 


27 


84 


55 


46 


15 


80 


41 


97 


25 


61 


± 





10 


24 


71 


± 





07 


24 


10 


± 





06 


1 


50 





90 





61 


4 


98 


-0 


13 


± 


0.41 


15 06 


26 


59 


55 


45 


05 


26 


42 


54 


22 


30 


± 





02 


21 


61 


± 





02 


21 


15 


± 





01 


1 


15 





69 





46 


2 


30 


-1 


36 


± 


0.18 


15 06 


28 


82 


55 


46 


50 


49 


42 


61 


25 


57 


± 





06 


24 


55 


± 





04 


24 


09 


± 





03 


1 


48 


1 


02 





46 


2 


42 


0.08 


± 


0.28 


15 06 


27 


64 


55 


45 


46 


80 


42 


66 


22 


18 


± 





03 


21 


42 


± 





02 


20 


91 


± 





02 


1 


27 





76 





52 


1 


99 


-0 


95 


± 


0.20 


15 06 


28 


45 


55 


46 


25 


36 


42 


86 


24 


10 


± 





03 


23 


47 


± 





03 


23 


02 


± 





03 


1 


08 





63 





45 


3 


40 


-1 


65 


± 


0.20 


15 06 


28 


71 


55 


46 


35 


43 


43 


84 


23 


98 


± 





03 


23 


26 


± 





02 


22 


80 


± 





02 


1 


18 





72 





46 


2 


97 


-1 


24 


± 


0.20 


15 06 


26 


69 


55 


44 


44 


60 


44 


97 


22 


17 


± 





02 


21 


50 


± 





02 


21 


05 


± 





01 


1 


13 





67 





45 


3 


19 


-1 


46 


± 


0.18 


15 06 


27 


85 


55 


45 


35 


12 


45 


04 


22 


90 


± 





03 


22 


12 


± 





02 


21 


65 


± 





02 


1 


25 





78 





47 


2 


21 


-0 


95 


± 


0.20 


15 06 


27 


46 


55 


45 


09 


04 


45 


43 


23 


60 


± 





03 


23 


02 


± 





02 


22 


58 


± 





02 


1 


02 





58 





44 


3 


08 


-1 


92 


± 


0.18 


15 06 


27 


97 


55 


45 


33 


02 


46 


01 


23 


94 


± 





04 


23 


00 


± 





03 


22 


49 


± 





03 


1 


45 





94 





51 


2 


02 


-0 


15 


± 


0.24 


15 06 


26 


83 


55 


44 


25 


46 


47 


61 


23 


47 


± 





03 


22 


68 


± 





02 


22 


20 


± 





02 


1 


27 





79 





48 


2 


93 


-0 


87 


± 


0.19 


15 06 


29 


35 


55 


46 


22 


64 


48 


30 


23 


47 


± 





03 


22 


74 


± 





02 


22 


27 


± 





02 


1 


20 





73 





47 


2 


71 


-1 


15 


± 


0.19 


15 06 


29 


46 


55 


46 


11 


11 


48 


43 


24 


20 


± 





04 


23 


42 


± 





03 


22 


93 


± 





03 


1 


27 





78 





49 


2 


96 


-0 


91 


± 


0.23 


15 06 


29 


03 


55 


45 


24 


94 


48 


95 


23 


54 


± 





03 


22 


87 


± 





03 


22 


40 


± 





02 


1 


14 





66 





48 


3 


24 


-1 


45 


± 


0.20 


15 06 


30 


53 


55 


46 


32 


97 


49 


34 


24 


70 


± 





04 


24 


01 


± 





03 


23 


54 


± 





03 


1 


17 





70 





47 


2 


71 


-1 


31 


± 


0.22 


15 06 


29 


94 


55 


45 


59 


62 


50 


20 


21 


84 


± 





03 


20 


96 


± 





02 


20 


44 


± 





02 


1 


40 





88 





52 


2 


29 


-0 


40 


± 


0.21 


15 06 


31 


60 




47 


08 


10 


50 


70 


25 


01 


± 





04 


24 


38 


± 





03 


24 


02 


± 





03 





99 





63 





36 


2 


60 


-1 


85 


± 


0.21 


15 06 


29 


61 


5 5 


45 


27 


53 


51 


80 


23 


03 


± 





03 


22 


22 


± 





03 


21 


65 


± 





02 


1 


38 





81 





57 


1 


99 


-0 


59 


± 


0.21 


15 06 


30 


61 


55 


46 


08 


54 


52 


50 


25 


13 


± 





07 


24 


46 


± 





05 


24 


01 


± 





05 


1 


12 





67 





45 


2 


12 


-1 


48 


± 


0.30 


15 06 


31 


71 


55 


46 


55 


76 


53 


60 


23 


42 


± 





02 


22 


87 


± 





02 


22 


42 


± 





02 


1 


00 





55 





44 


4 


49 


-2 


01 


± 


0.17 


15 06 


29 


24 


55 


44 


55 


99 


55 


55 


24 


74 


± 





04 


23 


92 


± 





03 


23 


42 


± 





03 


1 


31 





81 





50 


2 


02 


-0 


74 


± 


0.23 


15 06 


30 


34 


55 


45 


21 


76 


56 


10 


24 


01 


± 





05 


23 


17 


± 





04 


22 


67 


± 





04 


1 


34 





84 





50 


2 


02 


-0 


63 


± 


0.27 


15 06 


30 


41 


55 


45 


18 


26 


56 


60 


24 


23 


± 





06 


23 


46 


± 





04 


23 


02 


± 





04 


1 


21 





77 





43 


1 


96 


-1 


07 


± 


0.28 


15 06 


31 


41 


55 


45 


54 


54 


57 


10 


23 


67 


± 





05 


22 


80 


± 





04 


22 


31 


± 





04 


1 


35 





86 





49 


2 


33 


-0 


54 


± 


0.27 


15 06 


31 


75 


55 


45 


58 


62 


57 


75 


24 


17 


± 





05 


23 


56 


± 





04 


23 


08 


± 





04 


1 


09 





61 





48 


2 


77 


-1 


65 


± 


0.26 


15 06 


31 


15 


55 


45 


21 


98 


58 


45 


23 


98 


± 





04 


23 


29 


± 





03 


22 


82 


± 





03 


1 


16 





68 





47 


1 


65 


-1 


37 


± 


0.23 


15 06 


31 


43 


55 


45 


25 


88 


58 


45 


23 


28 


± 





04 


22 


42 


± 





03 


21 


93 


± 





03 


1 


34 





85 





49 


1 


97 


-0 


58 


± 


0.22 


15 06 


31 


95 


55 


45 


46 


99 


61 


60 


23 


93 


± 





04 


23 


13 


± 





04 


22 


62 


± 





03 


1 


31 





80 





52 


2 


74 


-0 


77 


± 


0.25 


15 06 


31 


04 


55 


44 


55 


54 


63 


15 


24 


18 


± 





04 


23 


42 


± 





03 


22 


95 


± 





02 


1 


23 





76 





47 


2 


44 


-1 


03 


± 


0.22 


15 06 


30 


85 


55 


44 


43 


28 


63 


30 


23 


64 


± 





03 


23 


00 


± 





02 


22 


55 


± 





02 


1 


10 





64 





46 


2 


63 


-1 


60 


± 


0.19 


15 06 


32 


99 


55 


46 


19 


41 


64 


35 


23 


78 


± 





03 


23 


10 


± 





02 


22 


63 


± 





02 


1 


15 





69 





46 


2 


40 


-1 


38 


± 


0.19 


15 06 


32 


54 


55 


45 


55 


30 


64 


60 


21 


58 


± 





02 


20 


86 


± 





02 


20 


38 


± 





01 


1 


20 





72 





48 


3 


02 


-1 


19 


± 


0.18 
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Table 1 — Continued 





R.A. 






Decl. 




Rgc 




B 








V 








R 






(B-R)o 


(B-V)o 


(V-R)o 






[Fc/H] 


15 


06 


30 


42 


55 


44 


11 


.48 


65.30 


26 


14 


± 


() 


.08 


25 


.39 


± 





.00 


24 


.97 


± 





.05 


1 


17 





.76 





41 


4 


.71 


-1 


20 


± 


0.34 


15 


06 


32 


90 


55 


46 


04 


02 


66.30 


25. 


52 


± 


o 


.07 


24 


.75 


± 


o 


.05 


24 


.24 


± 





.05 


1 


.28 





.77 





.51 


7 


.83 


-0 


.92 


± 


0.30 


15 


06 


31 


.47 


55 


44 


53 


.38 


66.55 


24 


66 


± 


() 


.04 


24 


.08 


± 


() 


03 


23 


.67 


± 





.03 


1 


.00 





.58 





42 


2 


.54 


-1 


.94 


± 


0.21 


15 


06 


32 


27 


55 


45 


20 


.77 


67.05 


24. 


18 


± 


o 


.05 


23 


41 


± 


o 


.04 


22 


.92 


± 


(1 


.03 


1 


.26 





.76 





.50 


5 


.66 


-0 


.97 


± 


0.25 


15 


06 


32 


.27 


55 


45 


20 


.77 


67.60 


24. 


18 


± 


o 


.05 


23 


.34 


± 


o 


.04 


22 


.84 


± 


{1 


.03 


1 


.33 





.83 





.50 


5 


.66 


-0 


05 


± 


0.25 


15 


06 


32 


.60 


55 


45 


10 


.58 


70.45 


23. 


82 


± 


o 


.04 


23 


00 


± 


o 


.03 


22 


.53 


± 


o 


.03 


1 


.29 





.76 





.52 


2 


.94 


-0 


.91 


± 


0.22 


15 


06 


32 


96 


55 


45 


21 


.88 


70.70 


24 


54 


± 


o 


.00 


23 


02 


± 


o 


.04 


23 


.10 


± 


o 


04 


1 


.44 





.92 





.52 


2 


.10 


-0 


22 


± 


0.28 


15 


06 


32 


.67 


55 


45 


04 


.23 


71.05 


23 


60 


± 


() 


.03 


22 


.71 


± 


o 


.02 


22 


.20 


± 


(1 


.02 


1 


.40 





.89 





.51 


4 


.69 


-0 


.38 


± 


0.21 


15 


06 


31 


51 


55 


44 


04 


.28 


71.85 


23 


32 


± 


o 


.02 


22 


.69 


± 


() 


.02 


22 


.22 


± 


o 


.02 


1 


.10 





.64 





.47 


2 


.63 


-1 


.59 


± 


0.18 


15 


06 


34 


.88 


55 


46 


15 


.02 


75.20 


24 


29 


± 


o 


.03 


23 


.65 


± 


o 


.03 


23 


.22 


± 


o 


.02 


1 


.07 





.64 





.44 


4 


.59 


-1 


05 


± 


0.20 


15 


06 


33 


.61 


55 


45 


12 


48 


75.60 


23 


83 


± 


o 


.04 


22 


90 


± 


o 


03 


22 


.44 


± 


o 


02 


1 


.39 





.93 





40 


2 


.08 


-0 


.30 


± 


0.23 


15 


06 


33 


.95 


55 


45 


26 


45 


78.15 


22. 


54 


± 





.03 


21 


.83 


± 





.02 


21 


.33 


± 





.02 


1 


.21 





.72 





49 


2 


.48 


-1 


.18 


± 


0.19 


15 


06 


36 


43 


55 


47 


12 


.95 


80.85 


23. 


87 


± 





.02 


23 


.14 


± 





.02 


22 


.67 


± 





.02 


1 


.21 





.73 





.47 


3 


.56 


-1 


.15 


± 


0.19 


15 


06 


33 


.72 


55 


44 


51 


40 


80.90 


26. 


17 


± 





.09 


25 


.38 


± 





.07 


24 


.84 


± 





.07 


1 


.33 





.79 





54 


3 


.56 


-0 


.76 


± 


0.40 


15 


06 


36 


77 


5 5 


47 


01 


61 


81.65 


23. 


90 


± 





.02 


23 


.16 


± 





.02 


22 


.65 


± 





.02 


1 


.25 





.74 





.51 


4 


.60 


-1 


.03 


± 


0.19 


15 


06 


36 


.03 


55 


46 


06 


.86 


82.20 


25. 


51 


± 





.05 


24 


.69 


± 





.04 


24 


17 


± 





04 


1 


.33 





.82 





52 


5 


11 


-0 


.68 


± 


0.26 


15 


06 


35 


.05 


55 


45 


16 


.93 


82.55 


22. 


90 


± 





.03 


21 


.95 


± 





.02 


21 


.42 


± 





.02 


1 


.48 





.95 





53 


1 


.92 


-0 


00 


± 


0.21 


15 


06 


35 


.14 


55 


45 


17 


.05 


84.35 


25. 


36 


± 





.08 


24 


.51 


± 





.06 


23 


.98 


± 





.05 


1 


.37 





.84 





53 


5 


.30 


-0 


.55 


± 


0.34 


15 


06 


34 


79 


55 


44 


57 


06 


84.70 


22. 


65 


± 





.03 


21 


.94 


± 





.02 


21 


.48 


± 





.02 


1 


.17 





.71 





40 


2 


.39 


-1 


20 


± 


0.19 


15 


06 


36 


.71 


55 


46 


24 


84 


85.40 


24. 


95 


± 





.04 


24 


.25 


± 





.03 


23 


.74 


± 





.03 


1 


.21 





.70 





.51 


1 


.90 


-1 


.22 


± 


0.22 


16 


06 


35 


.47 


55 


45 


06 


86 


85.70 


25. 


78 


± 





.09 


24 


.98 


± 





.07 


24 


.41 


± 





.06 


1 


.37 





.80 





.57 


3 


.53 


-0 


02 


± 


0.39 


15 


06 


36 


.17 


55 


45 


35 


06 


87.30 


25. 


05 


± 





.05 


24 


.10 


± 





.04 


23 


.40 


± 





.03 


1 


.65 





.94 





.71 


4 


22 


0.27 


± 


0.27 


15 


06 


36 


32 


55 


45 


34 


.78 


90.50 


24. 


18 


± 





.03 


23 


54 


± 





.03 


23 


.11 


± 





.03 


1 


.07 





.64 





.43 


2 


.06 


-1 


05 


± 


0.20 


15 


06 


36 


.06 


55 


45 


15 


93 


90.70 


24. 


11 


± 





.04 


23 


40 


± 





.03 


22 


.94 


± 





.03 


1 


.18 





.71 





.47 


2 


.25 


-1.27 


± 


0.22 


15 


06 


35 


.13 


55 


44 


24 


93 


91.85 


25. 


47 


± 





.05 


24 


.49 


± 





.03 


23 


.95 


± 





.03 


1 


.51 





.97 





.54 


4 


.98 


0.06 


± 


0.27 


15 


06 


35 


.30 


55 


44 


29 


.55 


92.75 


26. 


06 


± 





.08 


25 


22 


± 





.06 


24 


71 


± 





.05 


1 


.36 





.85 





.51 


3 


.11 


-0 


.58 


± 


0.34 


15 


06 


36 


.26 


55 


44 


43 


19 


93.05 


25. 


65 


± 





.06 


24 


.78 


± 





.04 


24 


.27 


± 





04 


1 


.38 





.87 





51 


6 


.44 


-0 


.47 


± 


0.30 


15 


06 


37 


.02 


55 


45 


12 


83 


93.40 


25. 


57 


± 





.07 


24 


.74 


± 





.05 


24 


.25 


± 





.05 


1 


.31 





.82 





.49 


3 


.97 


-0 


.72 


± 


0.32 


15 


06 


38 


.97 


55 


46 


40 


61 


95.15 


23. 


87 


± 





.02 


23 


.33 


± 





.02 


22 


.80 


± 





.02 


1 


.06 





.53 





.53 


4 


.48 


-1 


.95 


± 


0.18 


15 


06 


37 


.47 


55 


45 


23 


.02 


95.35 


23. 


89 


± 





.03 


23 


24 


± 





.02 


22 


.79 


± 





.02 


1 


.09 





.65 





.45 


2 


.37 


-1 


.58 


± 


0.20 


15 


06 


38 


.47 


55 


45 


50 


.92 


95.65 


26. 


22 


± 





.09 


25 


40 


± 





.07 


25 


.05 


± 





.06 


1 


.17 





.76 





41 


2 


.96 


-1 


.19 


± 


0.38 


15 


06 


36 


.56 


55 


43 


58 


.92 


95.85 


24. 


87 


± 





.04 


23 


.92 


± 





.03 


23 


.42 


± 





.02 


1 


.45 





.94 





.51 


4 


.51 


-0 


14 


± 


0.23 


15 


06 


39 


.13 


55 


45 


34 


77 


96.45 


26. 


31 


± 





.09 


25 


.54 


± 





.07 


25 


.03 


± 





.06 


1 


.28 





.78 





.50 


3 


.42 


-0 


.90 


± 


0.39 


15 


06 


39 


.01 


55 


45 


19 


.96 


97.55 


23. 


93 


± 





.03 


23 


.16 


± 





.02 


22 


.62 


± 





02 


1 


.31 





.77 





.54 


3 


.02 


-0 


.85 


± 


0.20 


15 


06 


39 


.55 


55 


45 


38 


10 
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